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Abstract —Analytical solution of Laplace equation using variable separation method, consists of two velocity 
potentials. However, only one component has been used. This research used both velocity potential equation 
components. With the potential equation, water wave surface equation and the related wave constants were 
formulated using kinematic free surface boundary condition and surface momentum equation. The characteristic of 
water wave surface that was produced was observed, both in deep water and shallow water. 
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I. INTRODUCTION 

The completion of Laplace equation using variable 
separation method in Dean (1991), produces two potential 
velocities, i.e. cosine and sine components. However, the 
application only exists in the cosine component in 
formulating various characters of water wave. This research 
does not discuss the use of the second component in Dean 
(1991), rather it studies the characteristic of water wave 
surface if the two velocity potential components are used 
simultaneously, i.e. water wave surface equation is 
formulated using a complete velocity potential. 

The formulation begins by formulating the final form of the 
two velocity potential components. Then, in each velocity 
potential, water wave surface and wave amplitude equations 
were formulated. Using wave amplitude equation and 
surface momentum equation, equations for wave number 
/{and wave constant G were formulated. 

It is obtained that the two velocity potential components 
have similar wave number k, wave constant G, but with 
different water wave surface equations. As the final water 
wave surface is the superposition or the sum of the two 
water wave surface equations. 

The characteristic of the water wave surface equation 
consists of maximum wave amplitude in a wave period, 
wave length, correlation between wave amplitude and wave 
height and other produced water wave surface profile was 
studied. 


II. TOTAL VELOCITY POTENTIAL 
EQUATION 

The form of Laplace equation solution (Dean(1991)), after 
periodic boundary condition was performed against time tis 
<pGc, z,t) = A cos kx(Ce kz + De~ kz ) sin at 

+B sin k x(C+ De~ kz ) sin a t 

.( 1 ) 

2n 

Wherecr = —is angular frequency, whereas Tis wave 
period.Constantfcwas obtained by performing lateral 
periodic boundary condition and wave number /c = y- was 

obtained where Lis wave length. Therefore,in (1) there is 
only one value of wave number Land one value of 
wavelength L. Based on the linear characteristic of Laplace 
equation, then(l) can be written as, 

<p(x, z, t ) = 0 A Gc, z, t) + <P B Gc, z, t) ....(2) 

Where, 

cp A (x, z,t ) = A cos k x(Ce kz + De~ kz ) sin a t (3) 

cp B (.x, z,t ) = B sin kx(.Ce kz + De~ kz ) sin o t (4) 

In (2), (3) and (4), the values of constants A, B, C and D 
should be determined. 

Equation (3) was performed at flat bottom (Dean (1991)) 
cp A Gc, z, f) = G a coshk (h + z) cos kxsino t (5) 

was obtained Similar procedure was performed in(4), 

< p B Gc, z, t) = G b coshk (h + z ) sin k xsin a t (6) 

was obtained 
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As has been mentioned, the two velocity potentials have 
similar wave number fc. There should have been one wave 
constant, i.e. G = G A = G B , but to ensure,a proof will be 
done in the following chapters. In the previous researches 
Hutahaean (2019 a, b) formulated equations for wave 
number fcand wave constant Gusing kinematic free surface 
boundary condition (KFSBC)and surface momentum 
equation. So is the case with this research, KFSBC equation 
and surface momentum equation will be used to formulate 
equations for wave number k and wave constant G. At the 
same time, this research is a improvementon the procedure 
of KFSBC integration against time t, in Hutahaean 
(2019a,b). 


III. THE FORMULATION OF WAVE NUMBER 
fcAND WAVE CONSTANT G A USING 0 A . 

3.1. Waterwave surface equation 

The first step in formulating equation for wave number 
fcand wave constant Gis the formulation of water wave 
surface equation to obtain wave amplitude equation. The 
formulation was performed using KFSBC. KFSBC equation 
using weighted total acceleration is (Hutahaean (2019 
a,b,c)). 


377 377 

y — = w - u — .(7) 

Whereyis weighting coefficient with the value of 2.87-3.14 
(Hutahaean (2019 c)). r)=r){x,t)\s water wave 

surfaceelevation against still water level (z = 0), is 
water particle velocity at horizontal-xdirection at the water 
surface (z = ij), whereasw^is the water particle velocity at 
vertical zdirection at the surface water. Using (5), equations 
of particles velocity at horizontal and vertical directions 
were obtained, i.e. 


u(x, z, t) 


3x 


= G A kcoshk (h + z)sinkx sinot...(8) 


w(x, z,t ) = — 


az 


= — G A ksinhk(,h + z)coskx sinot 


...(9) 

(8) and (9) were performed at z = tjand substituted to (7), 

y — = —G a ksinhk (h + rj) coskxsinot 
3 1 

—G A kcoshk 0i + rj) sinkxsinot ^...(10) 

Water wave surface equation was obtained by integrating 
(10) against time t. It’s visible that (10) is a non-linear 
function against time twhich is difficult to complete its 
integration. However, there are two arguments to make it 


simple, where the two arguments produce similar 
conclusion. 


The first argument is that the velocity potential equation 
was obtained using variable separation method, i.e. velocity 
potential d>that is regarded to have a form of &(x,z,t) = 
X(x)Z(z)T (t), whereVGc)is only a function of x, Z(z)is 
only a function of z andT(t)is only a function of timet. In 
this case Z(z) = coshkih + z). In relation with this, r )on 
the right side of the equation, both msinhk (h + r])or in 


coshkjh + rj) and—are not the function of time t, even 

3x 

though r) = r](x,t). Hence (10) can be written as, 

3rj 

at 

3m 


y~ = ~g a fc 


(coskxsinhk (h + rj) + sinkxcoshk (h + rj) si 
V 3x/ 


s inert 


.( 11 ) 


The second argument is that for a periodical function 
against timet , the element 

—G A k (coskxsinhk (h + rj) + 2 sinkxcoshk (.h + 

rj) — Jshould be a constant number against time t, which is 

strengthened with the formulation of velocity potential as a 
function of periodical time is just sinot. From the two 
arguments, the integration against time in (11) was 
completed only by integrating the stnertelement, obtained 

r, 

yo 

3m 


(coskxsinhk (h + rj) + sinkxcoshkjh + rj) —) 

v aw 


cosot 


•( 12 ) 


At the characteristics point, where in this research the 
characteristic point is a point where coskx = sinkx = 
cosot = sinot, (12) can be written as 


G A k 


1J Gc, f) = - 

yo 

(sinhk (h + rj) + coshkih + rj) coskxcosot .(13) 

The form coskx was selected because the first term of the 
elements in the parentheses is more dominant than the 
second element because of the presence of—element in the 
second term. It is defined a wave amplitude equation. 


A = (sinhk jh + rj) + coshk Or + rj) —).(14) 

yo V 3x/ 

Using (14) the waterwave surface equation becomes 
r)(.x,t) = Acoskxcosot 

Water wave surface equation was obtained at the velocity 
potential component^, i.e. 
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T) 0A (x,t) = Acoskxcosat .(15) 

^Voa , „ . , 

-= —kAsinkxcosat 

3x 

From (14), 

A„ 14 = — sinhk ( h + r] 0A ) .(16) 

•' 2 yo 

A v , 2A = ^ sinhk 0 1 + r ] 0A ) .(17) 

r] A (x, t) = ( A 711A coskx + A 71 2A sinkx)cosat (18) 

So,water wave surface equation consists of 4 (four) 
equations, i.e. (15), (16), (17) and (18), where wave 
amplitude in (15) is as input or known number. 

3.2. Equation for fcandG^. 

The next step is formulating equations for wave number 
/tand wave constant G A . The equation to calculate the two 
parameters can be obtained using (13) and surface 
momentum equation. (13) is differentiated against 
horizontal-* axis. 

a?7 G A k 2 


3* 


ya 


(sinhk (h + rf) + coshkfh + rf) sinkxcosat .(19) 

Bearing in mind that there are two variables that need to be 
calculated, then two equations are needed. As the second 
equation is surface momentum equation, where convective 
velocity is ignored. 

„ a 'f 


a u„ 


..( 20 ) 


3 Ujj 3 Uyj 

Substitute—where—was obtained from ( 8 ), and substitute 
at at 

—with(19), 

3% 


yaG A kcoshk(h + rf) sinkxcosat 

3rj 


(sinhk (h + g) + coshk 0 1 + g) —)st 
V 3x) 


sinkxcosat 


gG A k 2 

ya 

The equation is divided with G A kcoshk ih + 
g) sinkxcosat for sinkxcosat which is not the same with 
zero and remembering that in deep water tanhk (h + =1 

^-<4 -E) 

Substitute—with (19), 

3% 


i - 


G A k 2 

ya 


y 2 a 2 = gk 
(sinhk (h + rf) 


+ coshkQi + rf) —) sinkxcosat 
3x> / 


3rj\ 

— sir 

a x) ) 

Keeping in mind (14), 

y 2 a 2 = gk( 1 — kAsinkxcosat ) 


The left side of the equation is constant number, therefore 
the right side should be constant, maximum value of 
sinkxcosat = 1 is used 
y 2 a 2 = gk( 1 — kA) .(21) 

This equation is an equation to calculate wave number km 
the deep water. This equation has a maximum wave 
amplitude value and at the same time is a critical wave 
steepness for a wave period, i.e. in a large wave amplitude, 
(l — kA) = 0 can occur, or 
Amax = ^ .( 22 ) 

The calculation of (22) can be done if wavelength Lis 
already known. In the case that wavelength is not known, 
the equation for wave amplitude maximumcan be obtained 
by bearing in mind that (23) is a quadratic equation for 
wave number k, with a real root if the determinant is greater 
than zero. Wave amplitude maksimum was achieved at 
determinant valuedequal to zero, d = g 2 — 4gAy 2 a 2 = 0, 
so obtained 

..(23) 


By equating v4 maa .with (22) and (23), critical wavelength in 
a wave period was obtained, i.e. 

=^f~ 2 .(24) 

2 y z 


(22), (23) and (24) only apply for just one component, in 
this case is <t> A . 

As has been stated that from the two velocity potentials 
<t> A a.n&<t> B , there is only one value of wave number k, 
therefore it can be estimated that by using <t> B the form of 
wave number equation that is similar with ( 21 ) will be 
obtained. 

As an equation for G A , surface momentum equation (22) 
and water wave surface equation (15) were used and were 
performed at characteristic point. 

G a = -.(25) 


yacosk (k+“) 

IV. THE FORMULATION OF WAVE NUMBER 
fcAND WAVE CONSTANT C B USING0 c . 

4.1. Waterwave surface equation. 

Particle velocity equations to horizontal and vertical 
directions were formulated using , in ( 6 ). 


u(x, z,t) 


34 

3x 


= — G b kcoshk (h + z) coskxsinat 26) 


w(x, z, t) 


34 
az 
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= — G B ksinhk(.h + z)sinkxs inert... (27) 

Substitute (26) and (27) that was performed at z = r/to (7), 


and (36), where wave amplitude in (33) is as an input or 
known number. 


Y — = -G B ksinhk(h + rjjsinkxsinat 
3 1 

+ G B kcoshk(h + rjjcoskxsinat ^.(28) 

As has been performed in previous section, the right side of 
equation (28) can be written as, 

3 r] 

V 31~ 


— G B k [sinhk (h + rj)sinkx 

- coshkQi + rjtcoskx —) 
3x> 


5r l\ . 

— smut 

a x) 


am 

— coshkQi + rj)coskx —) cosat 
3x) 


...(29) 

Then, it was integrated against time t. 

r) (x, t) = 

G b k ( 

- (sinhk Qi + iysinkx 

ya V 

3 iV\ 

3x) 

....(30) 

At the characteristic point, (30) can be written as 
r) ( x, t) = 

G B k ( 3m 

- [sinhk Qi + ri) — coshkQi + n) — ) sinkxcosat 

ya \ 3x) 

.(31) 

Selected to use sinkx because the term in the parentheses on 
the right side is more dominant than the second term where 
there is Defined a wave amplitude equation, 

A = [sinhkQi + rj) — coshk ( h + r;) ....(32) 

Water wave surface equation becomes 

rj(x,t) = Asinkxcosat 

From (32) water wave surface equation was obtain from 
velocity potential component^ is 

t) ob ix, t ) = Asinkxcosat .(33) 

at?o b 


3x 


kAcoskxcosat 


IB 




ya 


sinhk Qi + ij 0A ) .(34) 


A r 2B = sinhk Oh + r] 0B ) —^— 

ya ,uc ax 

G b k 


,.(35) 


r] B (x, t ) = —-— ( A, 11B sinkx — A, l2A coskx) cosat 

.(36) 


Thus, water wave surface equation that was obtained with 
0 B also consists of 4 (four) equations, i.e. (33), (34), (35) 


4.2. Equation for fcandG B . 

(31)was differentiated against horisontal-x axis, 
a»7 G B k 2 / 

— =-1 sinhk (n + tj) 

3x ya V 

— coshkih + n) —) coskxcosat 

1 a x) 

.(37) 


Next, surface momentum equation was performed where 
convective velocity was ignored, 

.(38) 


du v a?? 

y — — ~q — 

r at ^ 3x 
a Ujj 

—was obtained from (26) whereas —from (37) and an 
at 3x 


equation was obtained, 

— yaG B kcoshkih + ri)coskxcosat = —g 

a?7 


G B k 2 


(sinhkQi +rj) — coshk Qi + rj) —) 
V 3x/ 


ya 

coskxcosat 


The two terms of the equation were divided with 
—G B kcoshkQi + rj)coskxcosat , and keeping in mind that in 
deep water tanhk Qi + rj) = 1 , 

3m 




Substitute(37) 


2 2 

y o 


ya 


gk 


[sinhk (h + rj) 

-coshkQi +r])—) 
3x) 


coskxcosatJ 


Bearing in mind (32) and by taking coskxcosat = 1, 
y 2 er 2 = gk(l - kA) ....(39) 

Equation for wave number k, i.e. (39) by usingd> B is the 
same as wave number equation formulated using <P A , i.e. 
(21), so it is proven that at potential velocity that is the 
superposition of two velocity potentials, each component 
has similar wave number . 

Furthermore, the equation for wave constant G B was 
obtained using surface momentum equation performed at 
characteristic point, obtained. 


gA 


yacosk (k+“) 


•(40) 


Compare to (25), the two velocity potentials have similar 
wave constant Gi.e.G^ = G B = G. 
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V. SUMMARY 

The description in Chapter III and Chapter IV shows that 
both velocity potentials have similar wave number equation 
and wave constant, so that both have similar wave number 
fcand wave constant G. Similarly, 


1. 


Wave number equation from (21) and (39) 
y 2 a 2 = gk — gAk 2 

Wave constant Gequation from (25) and (40), 
gA 


G = 


yacoshk (h + j ) 


2 . 


Water wave surface equation has a slightly 
different form. 

a. Water wave surface equation of 0 A 


rj 0A (x,t) = Acoskxcosat 
a?7 oa 


3x 


...(15) 

= —kAsinkxcosat 


A 


T),1A 


= —sinhkQi + g 0A ) .(16) 


ni,2A 


— sinhkih + .(17) 


a* 


r] A (x, t) = — {A 711A coskx + A r] 2A sinkx)cosat ....(18) 


b. Water wave surface equation ofd> B 

rj 0B ( x, t ) = Asinkxcosat .(33) 

3^°B , , 

■ = kAcoskxcosat 


5x 


_1 ?7,l B 


= — sinhk (h + g 0B ) .(34) 


ya 

Gk 


^11,2B 


sinhk Qi + t ] 0B ) 


3Vo B 


,.(35) 


V ^ (A v i B sinkx — A 71 2B coskx) cosat ....(36) 


c. Total water wave surface equation is 

h = h a+Vb 

VL RESULT OF MODEL 


6.1. The calculation of deep water depth 

Deep water depth was obtained using the criteria 


tanh k. 


o (h 0 + y) = 1. For a wave amplitude that is 

much smaller than deep water depth h 0 , k 0 h 0 (l + = 

k 0 h 0 = constant applies. tanhk 0 h 0 = lean be obtained at 
k 0 h 0 = a 0 7rwhere tanh (a 0 7r) = 1. Thus,fc 0 li 0 = a 0 7rand 

ho = f .( 41 ) 

Kq 

a 0 is a determined constant number, for example in CERC 


(1984) a 0 = 1 was used, so that — = 0.5 was obtained.hi 

this research, in addition to the criteria of tanh(a 0 n) = 
1, the value a 0 was also determined based on other reviews. 


Hutahaean (2019b) obtained that the larger the a 0 the larger 
the breaker depth and breaker height will be, so a 0 can be 
determined indiscriminately. 

For a large wave amplitude, tanhk (h + = 1, 

wherefe 0 ( h 0 + = an, revision on ar 0 was done 


againsta 0 , 

k 0 (h 0 -l 


..(42) 


Therefore the values of wave Ire u e n c y n a n d 4 0 p a ra me ter 
were absorbed in the value of a A . For the following 
calculation,a = a A was used. The value of <x 0 cannot be 
used too large, e.g. 2.25, wheretanh (2.257t) = 1, but it 
should take into consideration the characteristic of breaking 
that was produced, with the best value of a 0 = 1.6 — 1.9. 
Hutahaean (2019b)) obtained that with a 0 = 1.65, breaker 
depth that is in accordance with CERC (1984) was 
obtained. 


Table. 1: Wave characteristic in several wave periods 


T 

A 

n max 

£o 

h 0 

h o 

(sec.) 

(m) 

(m) 

(m) 

L 0 

6 

0,53 

6,9 

6,04 

0,88 

7 

0,72 

9,39 

8,22 

0,88 

8 

0,95 

12,26 

10,73 

0,88 

9 

1,2 

15,52 

13,58 

0,88 

10 

1,48 

19,16 

16,77 

0,88 

11 

1,79 

23,19 

20,29 

0,88 

12 

2,13 

27,6 

24,15 

0,88 

13 

2,5 

32,39 

28,34 

0,88 

14 

2,89 

37,56 

32,87 

0,88 

15 

3,32 

43,12 

37,73 

0,88 


Table (1) shows the result of the calculation of wave 
characteristic for several wave periods which includes deep 
water wave amplitude maximumA^*, deep water wave 
length L 0 anddeep water depth h 0 . The wave amplitude 
looks small but it will produce a large wave height, where 
the relation of wave height that is twice wave amplitude 
cannot be used. The calculation was done using the value of 
ar 0 = 1.75 the values of y = 2.05, where this value was 
obtained with the procedure in Hutahaean (2019 c,d), 
whereas as wave amplitude maximum, (23) was used. 
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6.2. Water wave surface profile 

The model was performed using wave period8 sec., wave 
amplitude 0.95 m . y = 2.05anda 0 = 1.75. The result of the 
model can be seen in Fig.l.a., Fig.l.b. andFig.l.c. 

Fig.l.a. shows that r] A andr] B (r] A water wave profile of 0 A , 
t) b water wave profile of 0 B ) have cnoidal profile, and both 
have similar profile size, i.e. wave crest elevation r) max = 
1.62 m, whereas wave trough elevation?j min = —0.66 m, 
therefore wave profile is asymmetric where T ) max is not the 
same as \r] min \. Wave height H = 1.62 + 0.66 = 2.28 m, 
Wilson parameter value (1963), T,m ^ x = 0.711, with this 
parameter value the wave profile belongs to cnoidal wave 
profile (Table (2)). 

Table.2: Wave type according to Wilson criteria (1963) 


Wave Type 

r l max 


H 

Airy waves 

< 0.505 

Stoke’s waves 

< 635 

Cnoidal waves 

^ r, r- ^Imax 

0.635 <-< 1 


H 

Solitary waves 

= 1 


thing that should be noticed is that there is a concavity in in 
the wave crest. A wave with a sharp wave crest can hardly 
be seen in a wave in the deep water, it always looks flat. 
Wave crest in Fig. l.b can be stated as flat, which is quite in 
accordance with the one in the nature. 


1.5 



x (m) 

- eta=etaA + etaB 

Fig.l.b. Wave profile rj = r] A + r] B in a wave period of 8 
sec., A = 0.95 m 



-1 

0 10 20 30 

x (m) 


- eta A -etaB 

Fig.l.a. Wave profile r] A andr] B in the wave period of 8 sec., 
A = 0.95 m 

The resultant wave, t] = t) a + t) b , (Fig.l.b.dan Fig. l.c.), 
obtainedj 7 max = 1.25 m, r/ min = —1.31 m, can be stated as 
symmetrical. Wave height H = 2.561 m, Wilson 
parameter ' ? "^ a * = 0.487, show that the wave has Airy’s 

wave profile type. The condition is very different from the 
ones previously known, i.e. Airy’s wave type can only be 
formed in a wave with a very small wave amplitude. One 


Fig.l.c. Wave profile in a wave period of 8 sec., A = 0.95 

in 

Furthermore, the calculation of water wave surface 
characteristic was performed in several wave periodswith 
wave amplitude maximum, i.e. in equation (37), A max = 
4 2 2 ■ The result of the calculation is presented in Table (3) 

and Table (4). 

Table.3: Water wave surface characteristics at wave 


amp litude maximum 


T 

(sec.) 

Vmin 

(m) 

Vmax 

(m) 

H 

(m) 

Vmax 

H 

8 

-1,31 

1,25 

2,561 

0,487 
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9 

-1,66 

1,58 

3,242 

0,487 

10 

-2,05 

1,95 

4,002 

0,487 

11 

-2,48 

2,36 

4,843 

0,487 

12 

-2,95 

2,81 

5,763 

0,487 

13 

-3,47 

3,3 

6,764 

0,487 

14 

-4,02 

3,82 

7,844 

0,487 

15 

-4,62 

4,39 

9,005 

0,487 


In the maximum wave amplitude, the Wilson parameter 
value ' lm ^ x = 0.487for all wave period shows that the wave 
belongs to to Airy wave, where?? ma;l . is quite close with 
I Vmin I- Furthermore in Table (4), the value of j = 
2.71 shows that an approach that the value of wave height H 
is twice the value of wave amplitude A cannot be 
determined or perfonned.Wave steepness ^ = 0.209, this 
condition exceeds the criteria of critical wave steepness 
from Michell (1893) with the value of^ = 0.142 

Table.4: The value of-^and wave steepness -Hn wave 
amplitude maximum. 


T 

A 

L 

H 

H 

(sec.) 

(m) 

(m) 

A 

L 

8 

0,945 

12,265 

2,71 

0,209 

9 

1,196 

15,522 

2,71 

0,209 

10 

1,477 

19,163 

2,71 

0,209 

11 

1,787 

23,188 

2,71 

0,209 

12 

2,127 

27,595 

2,71 

0,209 

13 

2,496 

32,386 

2,71 

0,209 

14 

2,894 

37,56 

2,71 

0,209 

15 

3,323 

43,118 

2,71 

0,209 


Note: wave height can be seen in Table (3) 

6.5. Water wave surface profile at breaker point 
To obtain water wave surface profile at breaker point, the 
values ofG, fcandTlare needed at breaker point. To obtain the 
value of the three wave parameters, shoaling and breaking 
analysis was performed. The shoaling and breaking model 
used in this research looks similar to the one in Hutahaean 
(2019 b), the model that was not discussed here. Bearing in 
mind that the two wave potentials have similar equations for 
wave amplitude, wave constant and wave number, then the 
shoaling and breaking model will also be similar to the 
model in Hutahaean (2019 b) that was formulated using & A . 


a. Water wave surface profile 

As an example of water wave surface profile at breaker 
point, a wave with wave period T =8 sec., wave amplitude 

A = 0.95 m and bottom slope— = —0.005was used. Water 

r dx 

wave surface profile at the breaker point is presented in Fig. 
2 a. and Fig.2b. 
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Fig. 2a. Water wave surface profile at breaker point, 
r], t] A andr] B 
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Fig. 2b. Water wave surface profile at breaker point ?? = 
V a + Vb- 


At the breaker point, wave profile is asymmetrical, where 
Pmax = 1-90 m, p min = —1.059 nr where wave height 
H = 2.958 nr whereas Wilson parameter 1,m ^ x = 

0.642 .With this parameter value, the type of the wave is 
cnoidal wave type. There are two phenomena that should be 
paid attention to, first, the occurrence of wave setup where 
in the deep water the wave profile is symmetrical, whereas 
in the shallow water the wave trough part is lifted. The next 
phenomenon is the separation of a wave from the two 
velocity potentials that were used, where there are two wave 
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crest. The presence of two adjacent waves also found in the 
coastal water. A more vivid example is tsunami wave on the 
coast or land, consist of two large main wave crests. 

In the profile of the breaking wave, it is also visible that 
there is a wave trough in fomt of wave crest. This also 
occurs in tsunami, where prior to the coming of the peak 
of the tsunami, the coastal water recedes first. 

b. Adjustment with breaker height index equation 

The result of breaker height model was calibrated against 
the average value of 5 (five) breaker height indexes. The 
adjustment was performed by multiplying wave constant 
Gwith 0.336. Whereas breaker depth h b was adjusted with 
breaker depth from SPM (1984), by changing the values of 
a 0 , wherea 0 = 1.76 resulted in a breaker height that fit 
with breaker depth from SPM (1984). The breaker height 
index equations used as comparators are breaker height 
index (BHI) equations from Komar and Gaughan (1972), 
Larson, M. and Kraus, N.C. (1989), Smith and Kraus 
(1990), Gourlay (1992) and Rattana Pitikonand Shibayama 
(2000), with the comparison result is presented in Table (5). 


Table.5: Comparison of breaker height model withBHl 


T 

(sec.) 

Wo 

(m) 

H b 

(m) 

(model) 

H b 

(m) 

(BHI) 

K 

(m) 

(model) 

h„ 

(m) 

(SPM) 

8 

2,602 

2,958 

2,955 

3,673 

3,697 

9 

3,293 

3,741 

3,74 

4,65 

4,678 

10 

4,065 

4,617 

4,617 

5,741 

5,776 

11 

4,919 

5,584 

5,586 

6,948 

6,989 

12 

5,854 

6,644 

6,648 

8,269 

8,317 

13 

6,87 

7,797 

7,803 

9,705 

9,761 

14 

7,968 

9,041 

9,049 

11,256 

11,321 

15 

9,146 

10,378 

10,388 

12,921 

12,996 


Note :BHI is the average value of breaker height from 5 
(five) breaker height index equations. 


As has been stated that the adjustment of breaker height was 
performed by multiplying wave constant Gat breaker 
pointwith a coefficient of 0.477. 

VIL CONCLUSION 

Both components of velocity potential equation as the 
solution of Laplace equation have similar wave number and 
wave constant, so that both can be performed as a unity to 
model water wave mechanics 


Water wave surface equation from each velocity potential 
component has different form, where the total of water 
wave surface equation is the sum of the two water wave 
surface equations. However, as has been stated that both 
have similar wave amplitude value and equation, wave 
number and wave constant G. Both produced similar wave 
profile. Therefore, both water wave surface equation are 
actually identical. Wave separation in the shallow water, 
also occurs in the nature, shows that the two velocity 
potentials should have been used, hi addition, water wave 
surface resultant have different wave height with each 
component of water wave surface. This also strengthens that 
the two velocity potential components of Laplace equation 
should have been used all simultaneously. 
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